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Abstract. Under conditions of clean flow we compute the leading term in the STF when 
the set of periods of the energy surface is discrete. Comparing to the case of non-degenerate 
periodic orbits, we obtain a supplementary term which is given in terms of the linearized 
flow. As particular cases, we give a STF for quadratic Hamiltonians and we obtain the Berry- 
Tabor formula for integrable systems. For conservative systems (i.e. systems with several first 
integrals), we give practical conditions to get a clean flow and interpret the leading term of 
the STF for a compact symmetry. We give several examples to illustrate our computation. 



1. Introduction 

Let H : M. 2n — > R be a classical smooth Hamiltonian. Under usual hypotheses (cf (|2.4p b one 
can define its Weyl quantization H and obtain a selfadjoint operator on L 2 (R n ). Let E GM. and 
e > 0. If H~ 1 ([E — e, E + s]) is compact then the spectrum of H is discrete in ]E — e, E + e[. 
To describe the spectrum in this interval, physicists introduce the local spectral density at E 
defined by the distribution 

(1-1) V E (h) = J2 s Mhh 

j>0 

where Xj(h) are the eigenvalues of H in \E — e,E + e[, and they study approximations of DeQi). 
On a mathematical point of view, it is more convenient to study the asymptotics when h — > + 
of a continuous fonction in E, 

(1.2) ~S S (/*) := ^Tr (^j(H)f 

where ip is a smooth compactly supported function in a neighbourhood of E (energy cut-off) and 
the Fourier transform of / is also with compact support. The function (|1.2p is the regularization 
of T>E{h) in a sense we should describe in section[5](cf Lemma |2~T|) . Asymptotics of Se(/i) when 
h — ► + is called the Semi-classical Trace Formula (STF in short). Modulo oscillatory terms, 
under some "clean flow conditions" (cf definition 12. 2p the expansion is a power series in h. Its 
remarkable property lies in the fact that coefficients of this series can be computed in terms of 
quantities describing periodic orbits of the classical Hamiltonian system associated to H in the 
energy shell He '■= {El = E}. The most celebrated version is the so called Gutzwiller formula for 
which one supposes that periodic orbits of He with period in a compact set are in finite number. 
The leading coefficient can be expressed as a sum over these periodic orbits which involves their 
primitive period, action, Maslov index and linearized Poincare map (cf (|3.2p ). 

This situation with only few periodic orbits appears for example for mixing flows, as the 
geodesic flow on a compact manifold of negative curvature. Nevertheless, in R™, it will never 
be satisfied when the Hamiltonian system owns some symmetries or first integrals. Indeed, if 
H is invariant by a one parameter symplectic group, then the image of a periodic orbit of He 
by an element of this group is also a periodic orbit of He of same period. This proscribes the 
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previous 'isolated' status of periodic orbits. Hence, the size of periodic orbit families grows with 
the symmetries. This is another motivation to investigate such systems, since the bigger periodic 
orbit manifolds arc in Eg, the more they participate to oscillations of Sb(/i) (see the leading 
power of h in Theorem I2.3|) . To our knowledge, only very few of such situations with symmetry 
were studied in this framework, which should be our concern in this paper. Our purpose is to 
give simple criterions to get a nice description of families of periodic orbits together with an 
asymptotic expansion of S_e(/i), and also to compute coefficients in such a way to interpret them 
as geometrical features of the classical motion as far as possible. 

Our article is structured as follows: in section[2l under quite general clean flow conditions, we 
perform a theoretical computation of the leading coefficient of the asymptotics in terms of the 
linearized flow and the structure of its algebraic eigenspaces (Theorem l2.4p . This is achieved by 
assuming that the set of concerned periods of Yje is discrete ('isochronous case'). Our calcula- 
tions are done pushing further a method based on coherent states due to Combescure, Ralston 
and Robert ([15j). which is to our opinion well fitted to make the linearized flow appear. In 
section |31 we apply this calculation to find back already known cases of STF, as the 'Weyl term' 
(zero period, cf (|3.ip ). the non-degenerate case ('isolated' periodic orbits, cf (|3.2p ). or the case of 
periodic flow on E^(cf (|3.3[) ) . We also give a STF for quadratic Hamiltonians (Proposition ^. 5p . 
In section^ we generalize the concept of non-degenerate periodic orbits to conservative systems 
(i.e. systems with several independant first integrals), via the notion of 'normal periodic orbits' 
already introduced in [34] and [38] . We show that it is a natural criterion to obtain a clean flow 
and a nice description of families of periodic orbits which arise in manifolds (Propositions 14.51 
and !4.8p . In section [5] we apply this last concept to the integrable case and interpret the leading 
coefficient in terms of action/angle variables. The STF involves the frequency of periodic tori 
and the Gaussian curvature of the energy shell in action coordinates (Theorem 15. 4p . This way, 
we recover the formula of physicists Berry and Tabor (cf [5J, 6 ) in a mathematical framework. 

Aknowledgements: We thank J. Bolte, B. Camus, A. Laptev, D. Robert and S. Vu Ngoc 
for stimulating discussions on the subject. 

2. Theoretical STF for isochronous periodic orbits 

2.1. Statement of the result. Let H : R 2 ™ — > R be a smooth Hamiltonian with associated 
dynamical system: 

(2.1) z t = JVH(zt), where J = 

We will use the notation z = (x, £) £ R" x R" as a variable in R 2 ™. The flow of H at time t G R 
with initial condition z € R 2 ™ is denoted by $t(z) = (qt,Pt)- The trajectory of z under this flow 
will be denoted by j z . If E £ R, then Y,e := {H = E} c R 2n . We define the monodromy matrix 
or linearized flow as 

(2.2) M g (t) := d g $ t (z) 

If A is in the spectrum of M z (t), wc define the algebraic eigenspace E\(z, t) by 

(2.3) E\ = E x {t 7 z) := £]ker(M,(i) - XId) k . 

k=l 

As we will be concerned with periodic orbits, the eigenvalue A = 1 should play a crucial role in 
the following. 

We make usual hypotheses of quantization on H which allow a nice functional calculus, 
namely, we suppose that there exists m > such that 

( <H(z)><C < H(z') >< z - z' > m , Vz,z'eR 2 ™. 

(2.4) <^ \d%H{z)\ < C a < H(z) >, Vz G R 2 " , Va € N 2 ™ . 
[ H has a lower bound on R 2 ™. 
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The Weyl quantization of H is denned as follows: for u € S(R"), 

(2.5) Opt{H)u(x) = {2*h)- n f f e^-^H(^Au{y)dydt 

In particular, under hypothesis (|2.4|) . Op^(H) is essentially sclfadjoint on S(R ra ) (see [17]), and 
we denote by D(H), H its selfadjoint extension on L 2 (R"). We fix E € R. We have in mind to 
study the spectrum of H near the energy E. For this purpose, we define the following regularized 
spectral density : 

(2.6) 5 E (h):=TrU(H)f 

where ip is smooth, compactly supported in a neighbourhood ]E — s, E + e[ of E (where e > 0) 
such that, if Eg := {H = E} C R 2 ™, then 

(2.7) H~ 1 ([E — e, E + e]) is compact and Eg has no critical point of H . 

ip(H) is an energy cut-off which is trace class by [27]. The function / is such that its Fourier 
transform / is compactly supported in R. The justification of definition (|2.6p is given by the 
following lemma whose proof is left to the reader: 

Lemma 2.1. Suppose moreover that < E\ < e and that ip = 1 on [E — E\,E + e{\. Then 
^T) E {h) = T) B {h) on (]£-£!,£+£![) and WD E (h))*fh = T;Se(h), where f T (t) := ±/(t/r). 

Thus, if ip is taken to be constant equal to 1 in a neighbourhood of E, and if L / = 1 (or 
/(0) = 1), then xSE(h) is simply a regularization of T>E(h) via a convolution by an approximate 
identity. The study of Sb(^) has proven to be an essential tool to investigate the spectrum of 
H. Well known applications are e.g. the asymptotics of the counting function of eigenvalues of 
H in an interval (see [37]), or the localization of the spectrum for integrable systems or periodic 
energy shells by Bohr-Sommerfcld conditions ([2D], [2], or [T2]V 

One condition to get asymptotics of B E (h) when h — > is to have a clean flow on supp/ xS^. 
Let us introduce the set 

G E ■= {(t,z) G supp/ x : $ t (z) = z). 

Definition 2.2. We say that the flow is clean on suppf xS^ if Se is a finite union of subman- 
ifolds o/R x R 2n and if at each point (T, z) of Ge, we have : 

{(r,a) G R x T Z S B : tJVH{z) + (M Z (T) - Id) a = 0} c T (t , z )Ge- 

When it holds, the last inclusion is actually automatically an equality. If the flow is clean on 
supp/ x Eb, then S_b(^) has an asymptotic expansion in powers of h when h — > + , possibly 
multiplied by oscillating terms of the form e^ a , a G R (see for example [3DJ, [7J, or [T5]). There 
are usual assumptions to ensure that the flow is clean. One can for example suppose that periodic 
orbits of supp/ x are non-degenerate. An other one is to suppose that the flow is periodic 
with same primitive period on Eg (see section [3]). We shall give later other explicit situations 
of clean flow, in particular for systems with several constants of motion (see section |4]). Let us 
denote by [Ge] the set of connected components of Ge- If (T, z) G Ge, then the classical action 

(2.8) A(T,z) := f Pt q t dt 

Jo 

is constant on each element of [Ge] (see [2J p. 167). The following theorem generalizes the 
results of Colin de Verdiere [TT], Chazarain |TU], and Duistermaat-Guillemin ([22]) on com- 
pact manifolds. Its proof goes back to the works of Guillemin-Uribe [25], Paul-Uribe [35] [36] , 
Brummelhuis-Uribe [7J, Meinrenken [2DJ or Dozias [TDJ, who gave mathematical sense to the 
heuristic first statement of physicists Gutzwiller [2B] and Balian-Bloch [?]. For more details on 
bibliography, we refer the reader to the nice survey |39] . 
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Theorem 2.3. Suppose that the flow is clean on suppf x Then we have when h — > + 

Sb(/»)= E {^h) 1 -^ L ei AY '^{ f f(t)d(t,z)da Y (t,z)+J2h j a j , Y ) +0(h°°). 

where day is the Riemannian measure on Y , and d(t, z) is a density to be specified on Y, we 
shall call the Duistermaat-Guillemin density ( cf [22] J . and Ay is the common action on Y. 

The Duistermaat-Guillemin density d(t, z) (in short DG-density) owns lots of the classical quan- 
tities characterizing periodic orbits of the Hamiltonian system (|2.1[) . Note that here it should be 
considered as a complex density, contrary to [52]. This is due to the fact that we should prove 
our results via coherent states and not using FIO. This method introduced in [H] skips the 
problems of caustics appearing in the usual WKB method, but requires complex phases, which 
explains the fact that our DG-density is also complex. Note that, to our point of view, the use of 
coherent states makes the computations more explicit in terms of the linearized flow M z (T) than 
the one given by FIO method (see Remark [T] of section |2~2"]) . As a consequence, our DG-density 
should include Maslov indices (see section [3]). Its computation is in general non-trivial and, to 
our knowledge, has only been achieved in some particular cases of importance (see section [3]). 
It usually appears as the determinant of a transversal Hessian coming from a stationary phase 
theorem. In [22] . Lemma 4.4, some computation of d is achieved in terms of half densities, when 
one has: 

(2.9) V(T, z) G G E , ker(M z (T) - Id) 2 n T Z ^ E = ker(M z (T) - Id) n T,H E . 

This hypothesis takes into account the cases of the period T — (Wcyl term), the case of non- 
degenerate periodic orbits or the case of periodical flow on Eg. However, when one considers 
systems with constants of motion in involution, this assumption is no more fulfilled as we shall 
see later (cf Lemma l2.5p . Our main motivation in this section is to broaden the computation to 
more general settings. Let us denote the set of periods by 

(2.10) L E := {T G supp/ : 3z G T, B , ^>t(z) = z). 
For T G L E , let 

Z T := {z£S £ : * T (z) = z}. 
As in |22j . we shall compute the DG-densities in the 'isochronous case\ i.e. when Le is finite. 

Theorem 2.4. Suppose that 

(2.11) L E is finite 

and that for all T in Le, Zt is a finite union of smooth manifolds such that Vz G %t 

(2.12) T Z Z T = ker(M z (T) - I) n T Z T, E and JVH(z) £ (M Z (T) - I)(T Z T, E ). 

Then the flow is clean on suppf x £# ; and we have modulo O(h +ao ) when h — > + : 
(2.13) 

5 E (h) = m J2 E (^hf-^ei^-^fiT)^- I f d(T,z)da Y (z)+'£h j a j , Y ] , 
TeL E Ye[z T ] ^ \ Y j>i J 

where [Zt] is the set of connected components of Zt, day is the Riemannian measure on Y , and 
A{T, Y) is the constant value of the action A on {T} x Y. 
Moreover if we suppose that 

(2.14) y(T,z)eG E , £i(T,z) = ker(A/ z (T)-M) 2 . 

then we have in addition that for all (T,z) in Ge 

(_l)«H fe+1 )det(u;o| E ) 

(2 15) d(T z) = v 

det(M z (T) - I) Wi \\n El ^H(z))\\ 2 det(U^J(M z (T)-I) lE5 )' 
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where u>o(., .) :=< J., . >r2™ is the usual symplectic form on R 2 ™, V\ = (JEi(T, z)) , and H El 
is the orthogonal projection on Ei(T,z). If £1 := ker(M z (T) — 7) n T Z H E , then k := dim£i ; 
£5 := £^ fl (Ei H T Z Y, E ), and ll£ 5 is the orthogonal projection on &s(T, z). 

Remarks 

• The determinant det (M Z (T) — I)\ v is non zero since, M Z (T) being a symplectic map, 
Vi(T,z) is the sum of all algebraic eigenspaces of M Z (T) corresponding to eigenvalues 
different from 1. 

• Note that there remains an ambiguity to obtain d(T,z) from d(T,z) 2 . This will give 
entire powers of e % ~ we should compute in particular cases, and this should involve 
Maslov indices of periodic orbits of Zt (see [25]). 

• We generalize the computations of [32J in the sense that (|2.9[) implies (|2.14p and JWH(z) ^ 
(M Z (T) — Id)(T z Y, E ). Under these supplementary assumptions in Theorem 12.31 we 
have £5 = {0} and 



(2.16) 



d(T,zf 



(-1) 



n+- 



det(w \ Bi ) 



det(M z (T)-I) ]vi \\U El (VH(z))\\ 



2.2. Proof of Theorem 12.41 We use the method employed in [TB] based on coherent states. 
We refer to [8] for the following reduction^ equations (3.15) and (3.17) of this reference). When 
h — > + , the trace S_e(/i) has the same behaviour in leading order than the following integral 



(2.17) 



I E (h) := 



(2nh) 
2V 



-d 



exp ( ~(p E (t,z) 



det I 



A t + iBt-i{C t +iD t ) 
2 



dtdz. 



where A t ,B t ,Ct, D t are the four n x n block matrices of 

A t B t 
C t D t 



M z (t) 



and ifE 



p2n 



is the complex phase given by 

•-PE = fl + ifl- 



(2.18) 



fi(t,z) := (E-H(z))t- \ I 



z)Jz s ds 



f 2 (t, z) := - < (I - Wt)(z t - z); (z t - z) > 



where Wt '■= 



W t -iW t 
-iW t -w t 

iB-t) • . Moreover, we have: 
(2.19) 

We send the reader to 



^ with ±(J + W t ) := (I - i l Mq)- 1 and M := (C t + iD t )(A t + 



\Wt\\„ m < 1. 



(Theorem 3.3) for the precise meaning of det* . We just recall that 
[det J] 2 = det. 

It is shown in [S] (Proposition 4.1) that if G E ■= {f' E = 0} n {%f E = 0} n (supp/ x R 2n ), then 



one has 
(2.20) 



G E = {(t, z) e supp/ x Yi E : $ t (z) = z}. 



indeed, assuming if x G K 2n is such that (M Z (T) - Id) 3 x = 0, then (M Z (T) - Id)x e T Z S B n 

ker(M z (T) - Id) 2 = T Z Y. E n ker(M z (T) - Id). Thus (M z (T) - Idfx = 0. This yields l2T4t . Moreover, if 
JVH(z) = (M Z (T) - Id)y, with y e T Z S B , then y e kcr (M Z (T) - Id) 2 nT z Z E = kei(M z (T) - Id) n T Z S E , thus 
J\7H(z) = 0, which is excluded. 

2 Take G = {Id} in this case. 
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(2.22) 



Moreover, an integration by part shows that 

(2.21) V(T,z)ee £ , <p E {T,z)=A(T,z). 

Remark 1. Let us compare with the usual WKB method using FIO. In this case, one boils down 
to the same kind of integral as in but with the real phase 

(fE{t, z) = S(t, z) — x.£ + tE, 

where S satisfies the Hamilton- Jacobi equation 

d t S{t,z) = -H(x,V x S(t,z)). 
S(0,z) = x.£. 

Actually our A(t, z) is just the constant value of S(t, z) — x.£ + tE on Ce (see e.g. |39j ). 

It appears that for the computations of the leading term of the asymptotics of S £?(/*), the phase 
tpE issued from coherent states seems to us more explicit in terms of the linearized flow M z (t) 
than the phase (f>E issued from FIO and defined implicitly via the equation V2.22]) . 

The conditions under which one can apply the generalized stationary phase theorem (see [15] 
Theorem 3.3) to integral (|2.17|) are exactly the clean flow conditions given by Definition ^. 21 (see 
[8] Proposition 4.3). Writing 

e E = |J m x Z T , 

it is easy to see that under hypotheses (|2.11|) and (|2.12|) . the flow is clean on (supp/) x 
Indeed in our case, if T e supp/ and &t(z) = z£ £e, then T(t,z)&e = {0} x T z Zt- Applying 
the stationary phase theorem yields modulo 0(h°°) when h — * + : 
(2.23) 

Q E (h) = m) E E (27rh) 1 ^ei A ^f(T)^- ( f d(T,z)da Y (z) + E ^ I . 
with 

(2.24) d(T, ,) := dot."* ( ^ z) ^A det ~i ( At + ^ " + U>t) 



where [Zt] is the set of connected components of Zt , day is the Riemannian measure on Y, 
and A(T, Y) is the constant value of the action A on {T} x Y . We have 

det ^ T,z)l »iT,. ) e E ^j = det f Hcss MT,z) + m TiTtZ) e E 

where ^t {t s) e E is the orthogonal projection on T(t. z )Ge and Hess lpe(T,z) is given by [8] 
Proposition 4.2 (take g = Id in our case). This implies that det ' J " (T, ' )t ' E 



/ | < (I - Wt)JVH(z); JVH(z) > 



det 



-±VH(z) 
V +\{ t M z {T)-I){I - Wt)JVH(z) 



-i\H(z) \ 
(*Af,(T) - /)(/ -Wt)JVH(z) 



2 



jr[JM z (T) + t (JM z (T))] 
+i(*M z (T) - J) (J - W T )(M Z (T) - I) + n £l / 



2 

where IT^ is the orthogonal projection on £i := T z Zt- Since M Z (T) is symplcctic, we have 
JM Z (T) +* (JM Z (T)) = CM Z (T) + I)J(M Z (T) - I). 

Set: 

(2.25) K := i(*M z (T) + I) J + i(*M,(T) -/)(/- W T ). 
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t X 1 J\IH{z) 


0VH{z) +* Xi (M z (T) — I) 


KJWH(z) 


K(M z (T)-I) + n El 



' t X 1 J\JH{z) 


i l WH(z) 4 


-* Xi(M z (T) — I) 


JVH(z) 


(M 2 (T) 




l -[(M z {T)-I) 


+ i(M z (T) 


+ i)4 



Then, the fourth block is equal to K(M Z (T) - I) + U £l . 

We recall that, classically, JVH(z) G kev(M z (T) - Id) and thus VH{z) G ker(*M z (T) - Id). 
Using this, we note that the third block is equal to KJVH(z). Let us set: 

(2.26) X 1 :=^{I-W T )JVH{z). 
We get: 

det , {T ' Z) E ' a = dot 

Now, in view of [8], Lemma 4.10 and our equation (|2.24[) . we have 
d(T,z)- 2 = (-l) n det 

Moreover^ 

(2.27) R- 1 = - 
We denote by 

(2.28) a :=< X u JVH{z) > . 

Note that a ^ since for a = (q,p), we have (|2.20p . < Wra,a >=< Wt(q — ip), (q — ip) >, and 
JVH(z) 7^ 0. We use the line operation L2 <— L2 — —JWH{z)L\, to get: 

(2.29) d(T, z)- 2 = (-l) n Q det(D). 
where 

D := (M Z (T) -I) + K^lle, - -JVH(z)[?VH(z) +* X X {M Z {T) - I)}. 

a 

Now, in view of l[2Tl"2"]) . we have X _1 n £l = |(M Z (T) + 1) JIL £l and therefore 

(2.30) £> = (M,(T) - J) - 1 -{M Z {T) + /)jn £l - i JV H {z^V H (z) +* X X {M,(T) - I)]. 

We are now going to compute the determinant of D in a basis of R 2n fitted to the linearized 
flow M Z (T). To do this we need a supplementary hypothesis describing E\(T, z). Let us suppose 
that we have 1(234]) . 

Lemma 2.5. [30] £ef M be a symplectic matrix and note E\ := X)fc>i kcr(Af — /<i) fe . Suppose 
that there exist two vector spaces V and W such that 

W GV CE 1 . 

V(x,y)£WxV, w o (x,y)=0. 
Then dim E\ > dim V + dim W . 

Proof. We recall that, if V\ := (JEi)^, then Vi is the sum of all algebraic vector spaces of M 
corresponding to eigenvalues different from one, and so we have Vi fl E\ = {0}. In particular, 
we get W l~l Vi = {0}, or equivalently, JW n = {0}. But, by hypothesis, we have JW C V 1 - 
and C V^. Thus 

dim F" 1 > dim E^ + dim W, 
which proves the result. □ 

Lemma 2.6. Let £2 = &2(T, z) := {a; G £1 : Vy G £ 1 , wo(a;, y) = 0}. VFe /lave 
Si = ker (M Z (T) - Id) 2 dimEi = dim£i + dim£ 2 . 



^Note that there is a slight misprint in [8] at this point, one minus sign is missing 
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Proof. Let us suppose that dim£i = dim£i + dim £2- Coming back to the proof of Lemma [2~5l 
with W = £2 and V = £1, we have for dimensional reasons 

(2.31) J£ 2 ®^ 1 L = £^. 

We recall that, by a symplectic argument, ker(* M Z (T) — Id) = Jkcr(A/ 2 (T) — Id). Thus we have 
(M Z (T) - Id)(E 1 ) C ( J Ex) 1 - H E x . Together with (pOTj) . this implies that 

(M Z (T) - Id)(E x ) c (£ 2 + (JEi) x ) r\Ei = £ 2 . 

Thus Ei = ker(M z (T) - Id) 2 . 

Reciprocally, let us suppose that Ei = ker(M 2 (T) — Id) 2 . By Lemma \2. 51 we have dimEi > 
dim£i + dim £2. Let us prove the inverse inequality. Let 3" be a, supplementary vector space of 
£1 in Ei. One has to show that dim^F < dim£ 2 . One of course have that (M Z (T) — Id)(3) + 
RJVH(z) C £ 2 . Thus 

(2.32) dim[(M 2 (T) - Id) (3) + RJVH{z)} < dim£ 2 . 

If kex(M 2 (T) - Id) C E E , then £j = kcr(A/ 2 (T) - Id) and dim(Af 2 (T) - Id)(3) = dimJ. In 
view of the last equation, we have dim? < dim £2. 

Now \{\ci{M z (T)-Id) is not included in E B , then cquivalcntly, JVH(z) g Im (M Z (T)-Id). 
Thus 

(2.33) dim[(A/ 2 (T) - Id) (3) + RJVH(z)} = dim(A/ 2 (T) - Id)(3) + 1. 

But, Y>e being an hypersurface of R 2 ™, there is at most one element of ker(M 2 (T) — Id) in 3. 
Thus dim J < dim (M Z (T) -Id){3) + 1. Together with ([2~52"1> and (jUSHJ), we get dim 3 < dim£ 2 , 
wheih proves our Lemma. □ 

We recall that we have K 2 ™ = Ei © Vi, where Vi — (JE1) 1 - and that Ei and Vi are invariant 
by M Z (T). We denote by 

(2.34) fc:=dim£i, r:=dim£ 2 . 

Lemma 2.7. There exists e x ^ in E x \ T Z E E such that (M Z (T) - Id)(ei) g RJVH(z). 

Proof. If ker(M 2 (T) - Id) C T Z Y, E , we have equivalently JWH(z) e Im (M Z (T) - Id). The 
hypothesis flTH?]) ensures that JVH(z) (M Z (T) - Id)(T z Y, E ), and we thus obtain e x . If 
ker(M 2 (T) - Id) is not included in T Z T, E , then we take e x in ker(M 2 (T) — Id) \ T Z T, E . □ 

We construct a basis of Ei as follows. Let 

• ai := JVH(z). 

• a2, . . • , a r be a basis of £2 

• a r+ \, . . . , as, be a supplementary basis of Span (ai, . . . , a r ) in £1. 

• ei be constructed as in Lemma 12.71 

• £2, ■ ■ ■ , £r be a supplementary basis of Span (ai, . . . , a^, £1) in Ei included in T Z T, E . 
For j = 1, . . . , k, we have 

D(a j ) = -^(M z (T)+I)J(a j ). 
Note that Vi C T Z Y, E since JVH(z) e E x . If v € Vi, then 

D(v) = (M Z {T) - Id)(v) --<Xi, (M Z (T) - Id)(v) > JVH(z) - Um z {T) + Id)m £l {v). 

a I 

If j = 2,...,k, then D(sj) = 

(M z (T)-Id){e 3 )--[< Xi,(M z (T)-Id)( Sj ) > +1 < VH{z),e 3 >]JVfr(«)-i(M,(T)+/d)jn El (e i )- 
a I 

D{ei) = — % - (M z (T) + Id).m £l (ei) --< VH(z),e x > JVH{z). 
2 a 

Now we can write the matrix of D in the preceedmg basis of M 2 ™. Note that since all D(ctj) 
(j = 1, . . . , k) span the image of |(M 2 (T) + I)JH £l , by column operations, we can remove its 
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contribution to all other columns. Then, since < VH(z),£i >^ 0, by substracting multiples of 
the column of D{e{) to other columns, we obtain a first line with only one non zero entry (the 
one of D(ei)). Note also that (M Z (T) — Id)(ej) <G £2, thus it can be expressed only in terms of 
qi, . . . , a r . Endly, we recall that (M Z (T) — Id)(v) <G V\. We use the following notations 

k r 

(2.35) 



~{M Z {T) + Id) J (a,) = J2 4<*i + E 



r • 



i=l 



where w 3 v G V\ ■ 



(2.36) 
(2.37) 



(M z (T)-Id)( Sj ) = ^^- 

i=l 

B := ((A}))2<i,j<r- 



det(D) = --< VH(z),si_ > det{M z (T) - Id)\ (-l) fcr (-l) fc dct(A) det(B) 



: = ((Wi))r+l<i,j<k+r, 

Then using all preceeding remarks, we obtain 

i 

a 

Taking the bracket of equation (|2.35|) with Ja p , for j,p = 1, . . . , k, wc obtain 

k r 

(2.39) i < aj,a p > = ^ w{ < cti, Ja p > + ^ w k+i < e ?; , Ja p > . 

i—r+l i—1 

Let X and Y be the k x k matrices defined by 

(2.40) Y := ((i < aj,a p >))i< P ,j<k 

< a r+ i, J a.\ > 



< a>k, Jai > 

< Oik 1 Jtt k > 



< El, JOL\ > 

<Si,Jotk> 



< s r , Joti > 

< S r , J(Xk > 



(2.41) X : = 

< a r +i, J®}, > 
Then (|2.39|) is equivalent to 

Y = XA. 

Using the definition of £2 and the fact that < ej, Jot\ >= for j = 2, . . . , r, we get that 

det(X) = (_i)''( fe -'')+i < ei,VH(z) > det(< a u Jaj >) r +i<i,j<k det(< e^Jctj >) 2 <i,j<r- 
We obtain, 



adet(.D) 



det(M z {T) - Id) Wi (-l) kr+1 i k+1 det(M)det(< a 3 ,a p >)i<j, P <k 
det(< Jai,Sj >) 2 <i,j<r{-l) r{k - r)+1 det(< on, Jaj >) r +i<i,j<k ' 



!J <r(-l) r(fe - r)+1 det(< on, Jctj >) r+ i 
Using ([235]) . this meanfl 

2 _ (-l)"dct(< Jcti,ej >) 2 <i,j< r det(< a;, Jo,- >) r +i<i,j<fc 



(2.42) 



d(T, z) 



i fe +! det(M z (T) - Jd)| Vi det(B) det(< a,, a 3 >)i<ij<k ' 

Taking the bracket of equation (|2.36p with —Je p , p = 2, . . . , fc, we obtain 

(2.43) det(< J{M Z (T) - Id)(e l ),e : j >) 2 <i tj < r = det(B) det(< Ja u e 3 >) 2 <i,j<v 

Lemma 2.8. // a 2 , ■ ■ ■ , ctk * s taken orthonormal and if e 2 ,...,e r is taken orthonormal and 
normal to E\, then 

\n. Ex {VH{z))\\ 



det{w \ Ei ) = det(< Ja t ,aj >) r +i<i,j<h 



■det(< ai,Jej >) 2 <i,j<r 



Proof. Compute the determinant in the orthonormal basis (ai, a 2l . . . , a k , £1, £ 2 , ■ ■ ■ , £ r ), where 
at := JVH(z)/\\JVH{z)\\ and := n £l (ViJ(z))/ ||n Bl (VH(z))||. □ 



Note that, since k + r = dimi?i, which is even, we have (— 1) 



r+k 
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Using (|2.43[) and (|2.42p , we obtain, under the assumptions of Lemma 
(-1)" det(w | El ) 



d(T,z 



2 



,-fe+l 



det(M z (T) - Id) Wi det(< J{M Z (T) - Id)(e ] ),e l >) 2 <i,j<r \\U El (\/H(z))\\ 
This ends the proof of Theorem 12.41 

3. Classical examples 



Let us briefly illustrate Theorem 12.41 for the most usual cases. In this section we only deal 
with examples such that (|2.9j) is satisfied. We let more general cases for sections 0] and EJ Thus, 
one can use formula (|2.16p . Of course, by a partition of unity, in the isochronous case, we can 
always suppose that supp/ PI £j E is reduced to a single point: 

3.1. The Weyl term. Suppose that supp/niL^ = {0}. Then locally G E = Y = {0} x S E . As 

M z (0) = Id, we have E 1 = R 2n , £i(0,z) = T Z Y, E , det(w | El ) = 1, U El (VH(z)) = VH(z), and 
Vi = {0}. Therefore 

d(0,z) 2 = 

||Vff(.)|| 2 

This case leads to asymptotics of the counting function of eigenvalues of if in a given interval of 
las/i goes to zero (see [37] ). Theoretically, one can also compute all the terms a^y of Theorem 
O using the asymptotics of Tt(ip(H)) for <p e C^°(]E - e,E + e[) described in (27J ("weak 
asymptotics"). It allows also to solve the ambiguity on the sign of d(0, z). We find in this case 



(3.1) 5 E (h)=TP(E)(^h)- n+l f(0)- ^. +0 {h- n+ % 

2n Jt, e l|Vi/|| 

where da^ E is the Ricmannian measure on T, E . 

3.2. Non degenerate periodic orbits. Let us recall that for a general Hamiltonian system of 

the type (|2.1|) . a periodic orbit (T,z) £ M* x R 2 ' 1 (i.e. such that $>t(z) — z). is non- degenerate 

if we have dimi?i(T, z) = 2. In a more geometrical way, 1 is not an eigenvalue of the energy 

reduced linearized Poincare map. This type of dynamics for periodic orbits arise in particular 

for the so-called 'mixing' systems (see the survey [IE])- One consequence is that {T} x 7 Z is 

an isolated connected component of 6^ (cf [3T] 'cylinder theorem' p. 136 Theorem 10), locally 

C E = {T} x 7 a . Since H El (VH(z)) is orthogonal to JS7H(z), we obtain an orthonormal basis 

r j-i j -1 , , ./ N \\n El (VH(z))\\ 2 

of Ei and easily compute det(w | El ) = ||vg(z)|| 2 

C— iW 

d(T,z) 2 - 



det(M z (T) - Id) Wi \\VH(zW 

Passing to d(T, z), one makes the Maslov index of the orbit appear, and after integration on the 
periodic orbit, this leads to a rigorous Gutzwillcr trace formula. Let us suppose that all periodic 
orbits of supp/ xE^ are non-degenerate. For sake of simplicity we suppose also that ^ supp/, 
in order to skip terms of zero period (cf (|3.1jl ). Then one has when h — > + 

(3 ' 2 ' SE( " ) = , ,L & 5 ^''^^w^m^m + ° , " ) ' 

where T* denotes the primitive period of 7, A(~f) is the classical action of 7 (.4(7) = J pdq), a 1 
is an integer called the Maslov index of 7 (see [Ml): and dP is the differential of the Poincare 
map of 7 restricted to its energy surface. 
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3.3. Periodic flow, (with unique primitive period). Suppose that T G K* is such that locally 
Ge = {T} x Tie- Differentiating < &r(z) = z on we obtain that M Z (T) = Id on T z E_e, and 
T z Yie C E\. As E\ is even dimensional (M Z (T) is symplectic), this yields E\ = R 2n . Moreover 
£1 = T z Yie- Thus 

d{T,z) 2 = -. 

\\VH{z)f 

This allows in particular to localize the spectrum in this setting, which yields Bohr-Sommerfeld 
conditions (see [20j). If all points of T*e have same primitive period T* then we have, without 
restriction on supp/, 

(3.3) SE(h) = m)(^hy n + 1 Y / hpT*)e^ A *e^±- f ^-+0(h- n + 2 ), 

where Ae is the classical action of any orbit of {Ae = f^pdq if 7 is an orbit of S_b), <je is 
their common Maslov index, and da^ E is the Riemannian measure on E^. 

3.4. Quadratic Hamiltonians. In this subsection, our aim is to show that already for the sim- 
ple class of quadratic Hamiltonians, periodic orbits can appear in families of arbitrary dimension 
and that this gives more 'exotic' STF. Here we restrict ourselves to the case where H is given 
by: 

(3.4) H w (x,0 ■■= \ m 2 + < Sx,x> Rn ), 

where S is a symmetric positive definite n x n real matrix. As S can be diagonalizcd in or- 
thonormal basis, it is straightforward to see that we can boil down to the case where S is 
diagonal, namely S = diag{w\, . . . , w^) where w = (u>i, . . . , w n ) € (IR^)™. Note that this system 
is "integrable almost everywhere" in the sense that the functions 

E(z):=±(tf+w 2 x 2 ), i=l,...,n 

are first integrals of the system which are in involution, but whose gradients are not linearly 
independant everywhere. Indeed, if z = (x,£) G R 2 ", whenever with have simultaneously Xj — 
£j = 0, the gradient of Fj vanishes at z and thus the gradients are colinear. We shall see that this 
situation precisely appears where the orbits arc periodic! As a consequence, in a neighbourhood 
of these points, the dynamics are not ruled by the Arnold-Liouville theorem (see [32]), since the 
level set of F := (Fi, . . . , F n ) are not necessarilly of dimension n . In the following, we investigate 
briefly several dynamics of the last type, corresponding to different diophantine assumptions on 
the WiS. We illustrate the preceeding remark showing that it is not necessary that periodic 
orbits arise in n-dimensional tori, and that they actually provide us with a various zoology of 
dynamics concerning families of periodic orbits. 

But before we recall some straightforward computations to solve (|2.ip very explicitly in our 
case. Here, the dynamical system is linear and we have 



cos(t v / 5) 


VS 1 sin(i v / S) 


-VS sin(tVS) 


COs(ty/S) 



M z (t) =:M(t) = 

So the flow is given by $t(z) = M(t)z, and the solutions of (|2.1|) by: 

( Xj(t) = cos(wjt)xj(0) + wj 1 smiwjfjtjiO). 
\ £j(t) = —Wj sm{wjt)Xj{Q) + cos{wjt)£,j{Q). 

Note that the set of periods of the dynamics is: 



(3.6) 
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So if T £ 9, then the set of all z for which T is a period is given by 
(3.7) A T := Re, +R e ;, 

where (ex, . . . , e n ; e' ± , . . . , e' n ) is the canonical basis of R 2 ". Note that, if z £ At, then VH(z) £ 
At- Thus At and E^ are transverse submanifolds in R 2 ". We can already notice that if 
q £ {1, . . . , n}, if we choose w\ = ■ ■ ■ = w q = 1 and all others Wj out of N, then dim(A27r flS^) = 
2<7 — l. This means that, within this family of quadratic Hamiltonians, we can find some C^; with 
isolated connected components of any dimension between 1 and 2n — 1 . Note also that after a 
short calculation , if z £ At, we have E\(T, z) = ker(M(T) — Id), and thus (|2.9p is satisfied. 

3.4.1. Different types of flows near periodic orbits. The proofs of the following lemmata are left 
to the reader. We start with the case owning the smallest number of periodic orbits. 

Lemma 3.1. The following assertions are equivalent: 

(1) Vi^j, f^Q. 

(2) All periodic orbits of H w are non-degenerate. 
Example: n = 2, w\ = 1, = \/2. 

Then we analyse the opposite case: 

Lemma 3.2. The following assertions are equivalent: 

(1) Vi,j,^£Q. ' 

(2) All orbits of H w are periodic. 

(3) 3T £ 9 \ {0} s.t. Vj, wjT £ 2ttZ. 

(4) 3z £ R 2ra s.t. rank(VFi(z), . . . , WF n (z)) = n, with z periodic point. 

(5) 3E > s. t. the flow is periodic on Tie- 

Example: n = 2, Wx = 1, W2 = 2, Note that here, periodic orbits can have different primitive 
periods. In our exemple, T\ = 2-7T, Ti — n are two different primitive periods. In this case we 
could have applied the Arnold-Liouvillc theorem. 

The very extrem is the case of the harmonic oscillator: 

Lemma 3.3. The following assertions are equivalent: 

(1) wi = •■■ = w n . 

(2) All orbits of H w are periodic with the same period. 

(3) 3E > s. t. on Tie oil orbits of H w are periodic with the same period.. 

Some cases in between are given in terms of symmetries, where periodic orbits should arise 
in manifolds of various dimensions. For a general Hamiltonian H we define the maximal group 
Gmax of symplectic symmetries by 

G = G max := {g £ Sp (n) s.t. H(gz) = H(z), Vz e R 2 "}- 

where Sp (n) is the group of symplectic matrices of size 2n. G max is obviously a subgroup 
of Sp(n). Let us suppose that it is a compact Lie group. The invariancc of H by the group 
classically implies that the flow $ f of H commutes with all elements of G max (these considerations 
are more detained in section l4~2|) . As a consequence, if (T, z) is a periodic point, then so is (T, gz). 
Thus, the whole manifold G{pf z ) is made of periodic orbits of same period and is included in 
Y>e where E := H(z). We now introduce a rough notion which ensures that these "tubes" of 
periodic orbits are "isolated" in I x E^. We should enlarge this notion in section [H If (T, z) is 
a periodic orbit of H, we say that (T, z) is " NDR in its energy surface" if we have: 

dim Ei (T, z) = dim(RJVH(z) + 3(z)) + 1. 

where S is the Lie algebra of G. This boils down to assuming that a periodic point in He becomes 
either a non-degenerate periodic orbit or a equilibrium point in Eg/G. When G max = {Id}, 
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we find back the notion of non-degenerate periodic orbit. The notation NDR stands for "Non- 
Degenerate in the Reduced space T,e/G" . In the case of the quadratic Hamiltonian H w , 

G m ax = Sp(n) n [{y/A^)- 1 0{2n)y/To\, where A = 

One can easily show that all periodic orbits are actually equilibrium points in the quotient. 
Thus, periodic orbit will appear by isolated 'tubes' of the form G(z) in 

Lemma 3.4. The following assertions are equivalent: 

(1) Si, 3, fjeQ^Wi =wj. 

(2) All periodic orbits of H w are NDR in their energy surface. 



o I) 



Example: d = 3, w\ — W2 — 1, W3 = \[2. 

3.4.2. STF for quadratic Hamiltonians. We use the notations of section [2l 

Proposition 3.5. For any w G (W.+ ) d , and any E in R* ; we have when h — > + 
(3.8) _ 

d 

Te U tt-z 

3=1 3 



2ir V2\l- cos(Tu)j)\ 1/2 JA T nz E \\^ H ( Z 



where At is given by [3.7\) , R(T) := card{j : WjT G 2-k'L}, dfTA r ns E *s the Riemannian 



measure 



on At nEfi, and ot is the common Maslov index of trajectories of At nSg. 
Proof. Le is finite since ? is discrete. Moreover, for all T in Le, we have seen that 
Z T = A T nS £; , VzeZ T , Ex(z,T) = A T = ker(M z (T) -7). 

Besides, H E and being transversal manifolds in R 2 ", Zt is a submanifold of R 2 ™ and Vz G 
Z T , T z Z t = kcr(M 2 (T) -/)nT z E B . We clearly have i[2~3]) . Hence, we can apply (|2~TS|) . 
An orthonormal basis of i?i is given by /? = Span {ej,e'j : Twj G 27rZ}. Since iuo(e£, e.,) = 
<5i.j, we obtain det(wo| J5 ) = 1. Performing a computation of eigenvalues of M Z (T), we obtain 
(e' ±m °'){i=i,...,»}- We get 



d{T,z) 



\VH{z)\Y II 2[l-cos(Tio^ 



Endly, the computation of A(T, Y) is also achieved using the expression of the flow. □ 

4. Conservative systems 

In this section we give practical means to obtain systems satisfying hypotheses of Theorem 
I2.4l for conservative systems, i.e. systems with several independent constants of motion near He- 

4.1. Normal periodic orbits. Let us note that if F : R 2n — > R is a first integral of the system 
(|2.ip , and if (To, zq) is a periodic orbit of He, then, differentiating 

F($ To (z)) = F(z), 

with respect to variable z at z , we obtain that VF(zo) G ker('M 2o (T ) — Id), and since M ZQ (T ) 
is symplectic, this means that 

(4.1) JVF(z ) G ker(M Zo (T ) - Id). 

This makes ker(M Zo (To) — Id) bigger for conservative systems. Moreover, E\ also grows dras- 
tically in this case as shows Lemma 12. 51 above all when first integrals are in involution. If 
our first integral F is independent of H at Zq, applying this lemma to M = M Zo (Tq), V = 
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Span (JVF(z ),JS7H(z )), W = RJV#(z ), we obtain that dim^i > 3. Thus, for conser- 
vative systems, we won't have non-degenerate periodic orbits. Another way to see this is to 
note that, if $f denotes the flow of F, then since it commutes with the one of H, all points 
of the image of* i ► $f (z ) are To-periodic of E^;. This is in contradiction with the fact that 
non-dcgcncratc periodic orbit E^i are isolated in R x E^. 

Our aim in this section is to broaden the notion of non-degenerate periodic orbit in order to 
be able to apply it to conservative systems. This notion should also ensure that the flow is clean 
on x supp/. For this, we make use of the concept of normal periodic orbits introduced in 
[31] and [HE]- Let E be a smooth manifold, and X a vector field on E with flow <f> t . Let (T, x) be 
a periodic orbit of X . Suppose that there exist first integrals F\, . . . ,Ff. in a neighbourhood of 
"/ x in E. Then, by using the same argument as in the beginning of this section, we always have 

k 

(M X (T) - Jd)(T x E) + RX(x) c f] keid x F r 

3=1 

Definition 4.1. We say that (T, x) is normal if there exist first integrals F\, . . . , F^ in a neigh- 
bourhood of "f x in E such that d x F\, . . . ,d x Fk are independent and that 

k 

(M X (T) - Id)(T x Y>) + RX(x) = f] ker d x Fj. 

3=1 

where M X {T) := d x $ T (x). 

A slight extension of the paper [3*4 to manifolds leads to the following 

Theorem 4.2. [34] Using the same notations, let (T,x) be a normal periodic orbit for X on 
E. Then there exists a neighbourhood U of 7^ in E, e > 0, and smooth submanifold 7 of 
}T - e,T + e[xU of dimension k + 1, such that (T',x') G]T - e,T + e[xU and ® T *(x') = x' if 
and only if (T', x') G CP. 

This is a kind of generalization of the 'cylinder theorem' (see [3T] or [32]) to normal periodic 
orbits. Note that in the case of the system (|2.1I) . using symplecticity of the linearized flow, 
(T, z) is normal periodic orbit in K 2 " if and only if there exists first integrals F\ , . . . , Fk 011 a 
neighbourhood of j z such that VFi(z), . . . , WFk(z) are independent and that 

(4.2) ker(M z (T) - Id) fl T Z E B = Span (JVF^z), . . . , JVF k {z)), 

and in general we always have 

Span (JVFi(z), . . . , JVF k (z)) C kcr(M 2 (T) -Id)nT z E B . 

We introduce a new definition close to the preceeding notion 

Definition 4.3. We say that {T,z) 6 R x E^ is "T,E-normal" if there exists first integrals 
Fi, . . . ,Ffr on a neighbourhood of ~f z in R 2n such that V-Fi(z), . . . , 'SZFk(z) are independent and 
that 

(M Z (T) - Id)(T z T, E ) +RJVH(z) = Span (VFi(z), . . . , Vi^z))-^ 2 " . 

Note that this notion is slightly different from saying that (T, z) is normal for the restriction 
of (|2.1[) to Eb since we suppose a little bit more assuming that our first integrals arc defined not 
only on E# but on an open set of R 2 ™. As a consequence, periodic orbits of the associated CP in 
E e will have the same period T (see Proposition I4.5[) . 

Lemma 4.4. (T, z) G Ge is Y,E-normal if and only if (T, z) is normal for the system 12.1)) in 
R 2 " and JVH(z) <£ (M Z (T) - Id)(T z Y) E )- 

Proof. If (T, z) g 1 x is E^-normal, then it is of course normal for the system (|2.ip in R 2 ™. 
Reasoning by contradiction, suppose that JVH(z) € (M Z (T) — Id)(T z T,E)- Let us denote by 
M the restriction of M Z (T) to T Z E^. Then rank(A/ — Id) = 2n — k, and by the rank theorem, 



SEMI-CLASSICAL TRACE FORMULA, ISOCHRONOUS CASE. APPLICATION TO CONSERVATIVE SYSTEMS 



dim(ker(M — Id)) = k — 1. This is in contradiction with the fact that ker(M — Id) already 
contains the fc-dimensional vector space generated by the vectors JV.Fi (2), . . . , JWFk(z). 

Reciprocally, suppose that (T, z) is normal for the system (|2.1|) in R 2 ™ and that JVH(z) ^ 
(M Z (T) - Id){T z Y> E ). Let W(z) := Span (VFx(z), . . .,VF k (z)). We have ker(M 2 (T) - Id) H 
T Z E E = JW{z). By the rank theorem, we get dim(M z (T) - Id){T z ll E ) = 2n — — 1, thus 
dim(M 2 (T) - Id){T z Y, E ) + JVH(z) = 2n - k. But since 

(M Z (T) - Id){T z j: E ) + RJWH{z) C (M Z {T) - M)M 2 " + RJWH(z) = W{z)^, 

we have equality for dimensional reasons. □ 

Proposition 4.5. Suppose that periodic orbits of suppf x T. E are all Y, E -normal. Then L E is 
finite and for all T £ L E , Zt has a finite number of connected components. These connected 
components are submanifolds of Y, E . Moreover the flow is clean on T, E x suppf, and for all 
T G L E , for all z £ Z T , JVH(z) <£ (M Z (T) - Id)(T z Z E ). 

Proof. We first prove a local lemma 

Lemma 4.6. If (T, z) £ G E and if Fx, . . . , are the k independent first integrals of Definition 
\4-S\ then 3e > 0, 3U E neighbourhood of j z in T, E , 3T E submanilfold ofU E such that 

e E n]T -e,T + e[xU E = {T} x T E . 

Moreover dimT E = k andVz' G T B , T Z >T E = JSpan {VF l {z'), . . .,WF k (z')). 

Proof. We can apply Theorcm l4.2l to the restriction of (|2.ip to T, E - Of course H doesn't count as a 
first integral. We obtain a neighbourhood U of z in Eg and e > such that Q E C]()T— e, T+e[xU) 
is a manifold of dimension k. If i £ {1, ... , k}, we denote by <& t i the flow of the Hamiltonian 
dynamical system generated by Fi. Note that since t ^ (T, $ t i (z)) is a smooth path in CbIHQT— 
e,T + e[) x [/), we have T (T , z) e E = {0} x Jiy( 2 ) where W(z) := Span (VFi(^), . . . , VF k (z)). 
Since the gradients of Fx, . . . , F k are still linearly independent close to z, for a smaller U, we can 
suppose that we have 

V(T', z') eG E n(}T-e,T + e[) x U), T (z ,. TI) £ E = {0} x JW(z'). 

Using a chart of G E D (]T — e,T + e[) x U) defined on a connected domain, it is then easy to 
show that, schrinking ]T — e,T + e[) x U, we have that C E n (Ux]T — e,T + e[) is of the form 
{T} x T E . □ 

Then, using the compacity of supp/ x S^;, we obtain easily that H E is finite (reason by 
contradiction). The set Zt satisfies the property Va; G Zt, 3k x £ N, 3U X open set of S^; such 
that Zt fl U x is a submanifold of dimension fc x of S^;. From this it follows easily that each 
connected component of Zt is a submanifold of T, E . Now to show that Zt has only a finite 
number of connected components, let us suppose by contradiction that we have a sequence 
(C n ) rl6 N of disjoint connected components of Zt- Pick a z n on each C n . Then by compacity 
of Ti E and thus of Zt, by extracting a subsequence, one can always suppose that z n converges 
to z* G Zt- Applying the preceedmg property to z* , one can always suppose that Zt H U z * 
is connected. Thus for n great enough, z„ £ U z * , and C n is the connected component of Zt 
containing z* . This contradicts the fact that C n are all disjoints. Endly, Lemma T4.6I garanties 
us that V(T, z) £ C E , we have £i(T, z) = T Z T E . Thus the flow is clean on T, E x supp/. □ 

We see that under normal periodic orbits assumptions, the hypotheses of the first part of 
Theorem 12.41 are fulfilled. Note that a non-degenerate periodic orbit (z,T) is E^-normal. We 
have even a more general statement: 

Lemma 4.7. Let (T,z) be periodic orbit ofY, E such that there exist independent first integrals 
Fx, . . . ,F k is a neighbourhood of 7 2 in K 2 " and such that dim Ex(T, z) = k + 1. Then (T, z) is a 
Y, E -normal periodic orbit. 

An example of such a situation is given by Lemma 13.41 
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Proof. First note that, without any normality assumption, E\ is never included in Eg. Indeed, 
Ei C E B <=> JVH(z) G (J-Ei) 1 - = Vi, and we have seen that EiHVi = {0}. Now, under our 
assumptions, we have thus dim(i?i (~l E#) = k and 

Span (JVJiO), . . . , JVF fe (z)) C kcr(M z (T) - Id) H T Z E B c£iH T Z E £ . 

For dimensional reasons, we have equality everywhere and thus (T, z) is normal in R 2n . Moreover, 
if a G T Z E B is such that JVH(z) = (M Z (T)-Id)a, then a G £inE B C hsr{M z {T)-Id)r\T z Y, E , 
thus JVH(z) = 0, which is excluded. The lemma is proved using Lemma [4~4l □ 

Note that, even if the assumption 'dim£' = k + V seems a more straightforward generalization 
of the concept of non-degenerate periodic orbit, in view of Lemma 12.51 it will actually not be 
fulfilled each time that the system has first integrals different of H in involution, i.e. satisfying 
{Fi, Fj} = 0, where {, } is the Poisson bracket. 

To end this subsection, we give an example of a periodic orbit which in not normal. Let us 
consider the quadratic Hamiltonian of section [3] in dimension n = 2, with wi = 1, and W2 = 2. 
In this case, we have two periods, tt and 2tt, for which A T = Re 2 + Re 2 and = R 4 . If we 
consider the period T = 2ir with z G A,r, we have $>t(z) = z, and, for any first integral F around 
j x , we have JVF(z) eA, = Ke 2 + Re 2 . But we have £i(T,z) = T Z E£;. Thus (JOJ) cannot be 
fulfilled for any system of first integrals. 



4.2. Compact symmetries. A nice way to obtain first integrals of Hamiltonian dynamical 
systems is due to E. Noether: one can look for the symmetries of the Hamiltonian H : M. 2n — > R. 
We suppose that there is a compact Lie group G C Sp (njf| such that 

(4.3) V.g G G, H{gz) = H(z). 

Actually, we should be concerned by groups of symplectic diffcomorphisms and not only linear 
groups. We do this for sake of simplicity. Note that since a group of isometries is affine, with our 
compactness assumption, we arc not far from reaching all symmetries of the problem. Moreover, 
for linear groups, Noether first integrals have an explicit expression: if A G S the Lie algcgra of 
G, then 

(4.4) F A (z) :=< JAz,z > R 2„, 

is a first integral of the system (|2.ip . Indeed, by derivating with respect to variable t G R at t = 
the identity H(e tA z) = H{z), we obtain {H, Fa} = since J A is symmetric. Now, if g G G, 
differentiating the identity (|4.3p with respect to variable z, we obtain, using the symplccticity of 
G that 

Vz G R 2n , gJVH(z) = JVH{gz). 

As a consequence, the flow of H commutes with all elements of G. Now, if (T, z) G R x Y E is 
a periodic orbit, then for all g in G, so is (T,gz). We obtain a manifold {T} x G(7 Z ) included 
in Ge- Note that when our group is not finite, this proscribes non-degenerate periodic orbits. 
We have in mind to ensure that {T} x G(7 Z ) is a connected component of Ce and that we get 
a clean flow at (T, z). In view of definition 12.21 a necessary condition for that is: 



(4.5) 



ker(M z (T) - Id) nT z Z E = RJVH{z) + Sz. 
JVH{z) i {M Z {T) - Id)T z E E . 



Morally speaking, this means that G has to be big enough so that WJWH(z) + Qz would 'fill' 
ker (M Z (T) — Id) n T z Es, or from an other point of view, that ker(M z (T) — Id) must not be too 
big. This condition is actually enough to ensure a clean flow as states the following 



'Sp(n) is the group of symplectic linear maps of R 2n . 
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Proposition 4.8. Suppose that H admits a compact connected Lie group of symmetry G such 
that Ge satisfies J^.5|). Then there exist periodic points z\,...,z r of Tie such that 



(4-6) e £ = |J |J «jxG( 72i ), 



i=l n£Z: 

nT* . £suppf 



where T*. is the primitive period of Zi. Moreover the flow is clean on suppf xEg, and if {2.1$ 
is fulfilled on suppf x Tie, then, assuming for simplicity that ^ suppf , one has when h — > + 



i—1 n£ 



nT*, Gsuppf 



^\&et{M Zt {TlY - Id) Wi \ 



(47) x f |det(«;o| Ei )^ ^^^^+1, 

■W,) ||n Bl (vir(z))|| ^|dct(n £5 j(M 2 (T)-/d)| £s )| 

where ki = dimG(7 2i ), Aj := J pdg, c ^ f7 G(7- )(z) * s ^ e Riemannian measure on G(~f Zi ] 
Z, and other terms are defined in Theorem \2.3\ 



Proof. From Proposition 14.51 and Lemma |4.6| we get that He is finite, the flow is clean and for 
all (T, z) in Ge, dim(ZT) = dim(G(7 2 )). As G(7 2 ) is a compact submanifold of same dimension 
as Zt, it is one of is connected components. This yields (|4.6[) . For the rest of the theorem, 
one applies Theorem 12.41 noting that if n £ Z, if z £ G(j Zi ), then the matrices M z (nT Zi ) and 
M Zi (T Zi ) n are conjugated. Indeed, if g £ G and t £ K, we have 

M gz (T) - gM z (T)g-\ M z (T) n = M z (nT) , M» tW (T) = A/ 2 (i)M z (T)M 2 (t)- 1 . 

This is obtained differentiating with respect to the z variable the identities <f>t+ s (z) = $ s (<I> t (z)) 
and g$>t{z) = <&t{gz). □ 



Such a situation is of course fulfilled in Lemma l3.4l An example where it is not fulfilled is the 
same as in the previous subsection. If we compare formula (|4.7|) to the one of the non-degenerate 
case (see p.2p ). we see that we obtain similar quantities, unless for the term T* which is here 
replaced by our integral on G(7 2i ). Of course, the ideal result would be to interpret this integral 
in terms of geometrical quantities describing G(7 2i ). A heuristical interpretation can be found in 
the work of physicists Creagh and Littlejohn ( [U] , [T7] ) . So far we have not been able to achieve 
this in a rigorous general mathematical framework for the moment. Nevertheless, we do so in 
the next section for a class of integrable systems. 



5. Integrable case: the formula of Berry- Tabor 

In this section we give a trace formula for a certain class of integrable systems in terms of 
action-angle variables. In this context, our result gives actually less informations on the spectrum 
than the Bohr-Sommerfeld conditions (see for example [13] [20] [2]) which allow to localize the 
spectrum, whereas we are dealing with spectral statistics. Nevertheless, our goal is to highlight 
the geometrical interpretation of the STF, in order to be able to do so in more general contexts, 
where less is known about the spectrum. Another paper by Charbonncl-Popov [5] may seem very 
close to what we are presenting here. However, the object of study is not the same since they 
are considering trace formulae defined with Schwartz functions in R fc of k commuting operators, 
which can not be reduced to a Schwartz function of only one of them. We should here recover 
results of physicists Berry and Tabor ([S][5]) in a rigorous way. To our knowledge, these results 
seems to be strangely absent from the mathematical literarure. 
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5.1. Assumptions and notations on the integrable system. Let H : R 2n — ► R be a 

smooth Hamiltonian, E £~R. Suppose that there exists £q > such that — £q, E + Eq)) is 

compact without critical points of H and that the system is Liouville integrable (or integrable) 
on Ue '■= H~ 1 (\E — £o,E + So[), i.e. assume that there exist n smooth constants of motion, 
Fj : R 2 ' 1 — > R, j = 1, . . . , n, such that we have: 



3 • 

(1) = 0, .'/../ ! //. on[/ £ . 

(2) VFi(z), . . . , VF„(z) are independent for all z in Ue- 



(5.1) 



where {, } is the Poisson bracket. We can always assume that F\ = H . We denote the flow of 

F ■ 

Fj by $ f 3 . These different flows commute since the integrals are in involution. Let 

F:= (F 1 ,...,F n ):U E ^R n . 

and if c £ R", T c := F _1 ({c}). Due to the independence of the first integrals, each T c is a 
compact submanifold of dimension n of R 2 ™. We will suppose that 

(5.2) For all c in R™, T c is connected. 

According to Arnold [3], each T c is diffcomorphic to a torus and it is covered by the joint flow 
of the Fj. We obtain also local action-angle coordinates near each T c (see [32]). Nevertheless, 
we will suppose that we have global action-angle coordinates on Ue- To this purpose we make 
the following hypothesis 

(5.3) ¥(Ue) is simply connected in R". 

According to Duistermaat [2T], we have global action-angle coordinates on Ue, i-e. there exists 
a symplectic diffcomorphism ip : Ue — ► T" x D, where T n := R ra /(27rZ)™, D is an open set in R™ 
such that, if we write tjj = (ifti j ^2) 

• \fz,z'eU E , ¥(z)=¥(z') ^ Mz) = M* r )- 

• Uip t = (Ot, It) ■= t/jo<£> t o tp-\ then for all (6, 1) £ T n x D we have 

j 6 t (6J) = VH{I). 
l ' j \i t (6J)=0. 

where H : D — > R is the smooth function defined by 

whatever in , since H is constant on the tori. Thus ip transform symplectically the motion 
into translations on tori. In other words, the flow ipt is given by: 

(5.5) <p t {6,I) = (6 + [tVH(I)],I) 

where [.] : R" — * T" is the canonical projection. One usually denotes by w : D — > R" the 
frequency function defined by 

(5.6) w(I) := V-ff(I). 

We suppose that our integrable system is non- degenerate, i.e. 

(5.7) V/ £ D, det(u/(I))) 7^ 0. 

This assumption is 'gencrically' satisfied and it is also often used for perturbation of integrable 
systems (see [3]). Let us now focus on periodic orbits of the system. Of course if one point of 
a torus is periodic, then all are. We will talk about 'periodic tori'. From (|5.5[) we deduce 
that for z in U E , T e R, we have § T {z) = z if and only if Tw(I) £ (2irZ) n , where ip 2 (z) =: I. 
Using this remark and assumption (|5.7|) , one can apply implicit function theorem to the function 
/ : R x D -> R™ defined by f(T, I) = Tw(I), to prove that, if (T , 1 ) is a periodic torus, then 
the set of periodic tori around (To,Io) in R x D is locally given by the graph of a function 
3 :]T - e, T + e[-> D such that 

(5.8) VTe]T -e,T + e[, T W (3(T)) = T u.(/ )- 
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We introduce 

(5.9) t E := {H = E}c R n . 

Note that since our dynamical system has no equilibrium point on E#, we have that VI G S^, 
VH(I) 7^ 0, and thus E# is a hypersurface of R". To enter the context of Theorem l2.4( we need 
to have a constant period on S^. Thus we would like the image of 3 to be transversal to E^ at 
Iq in R™, for which it is enough to have < 3'(Iq), VH(Iq) >^ 0. But derivating equation l|5.8p . 
we obtain that J' (To) = —^w'(Io)^ 1 {w(Iq)). This leads us to make the following hypothesis 

(5.10) V(T, J) G supp/ x T, E such that Tw{I) G (27rZ)™, < iy'(J)- 1 (w(7)), u;(J) >^ 0. 
We have proved the following 

Lemma 5.1. Make assumptions \5. 7j ) and \5.1U\) . Let (Tq^zq) G suppf xT,e such that <&t (zo) = 
zq. Then it exists e > suc/i that if (T, 7) G suppf x Eg, if |T — Tq\ + |7 — Iq\ < e and if (T, I) 
is a periodic torus then I = Iq and T — Tq. As a corollary, there exists only a finite number of 
periodic tori o/E^; with period in suppf. 

The STF formula should involve the function K : Es — > R of Gaussian curvuture of E^. We 
recall (sec 33 p.3B-21 5. Corollary) that if I G E#, if ip is a paramctrization of E^ defined on a 
neighbourhood of zero in R n_1 such that ^(0) = I, then K (I) is defined by 

detf< imp. ,#g-(0) >1 

(5.11) /,(/)- 1 |V*ttl'W / 1<<J<„-1 



det(<%^(0),c^(0) >) x 



<i,j/<n— 1 



The following lemma shows that hypotheses (|5.7p and (|5 . 1 0[) are actually mandatory if we want 
the curvature to be non-zero. 

Lemma 5.2. For all I G E^;, 

(5.12) detKoo) < u/urvoow^o >= K^^x-ir'Voor -1 - 

Proof. Let <p be a chart of E^ at I as in the definition (|5.1ip . Differentiating at £ = the 
relation 

(5.13) <«MO)>fl£ 4 v(0)>=0, 

we obtain that 
(5.14) 

det(< <//(!)% 93(0), %p(0) >)i<ij<„-i = (-l)"- 1 dct(< w(J),%d £i ^(0) >)i<ij<„_i. 

Let a := w'(I)~ 1 (w(I)). First suppose that a G T/E^. Then there exist real numbers 
fxi, . . . jJL n —i such that 



i.e. 



n-l 

= 5>X(iX0&¥>(o)). 

3=1 

Taking the bracket of the last equation with 9^(^(0), we obtain that K(I) = via (|5. 14[> . and 
equality (|5 . 1 2[) holds since in this case < w'(I)~ 1 (w(I)),w(I) >= 0. 

Suppose now that a ^ E#. Then e := (a, d^ 1 (p(0) 7 . . . , d^ n _ 1 ip(0)) is a basis of R". In view of 
equation (|5.13[) . so is e' := (w(I),d^ 1 ip(0), . . . , d^ n _ 1 Lp(0j). We write 

(5.15) det(u/(J)) = det(u/(J)) e , e ' det(/d)^, e . 
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Wc easily get det(Id) E > ;E = <lu >(/)-i"(^(jj) w (i)> ' Moreover, if a^j are real numbers such that 

n-l 

r=l 

then det(u/(/)) £i£ ' = dct(a;.j)i<ij<„_i. Tacking the bracket of the last equation with d^ip(0), 
we obtain that 

dct(< w'{I)d^{0),9^(°) >)=det{w'(I)) s , e , det(< %^(0),%<^(0) >)i< w <„-i. 
Together with equations (|5.14[) . and (|5. 15[) . this proves the Lemma. □ 

5.2. The Berry- Tabor formula. Our goal here is to use Theorem 12.41 under precccding hy- 
potheses to give a STF in the intcgrable case involving geometric quantities. The following 
lemma shows that assumptions of this theorem are satisfied. 

Lemma 5.3. Make hypotheses i5.1)) , I5.ty) , i5.3\) . Let z <G Sb, and denote by (8,1) = ip(z). 
Suppose that <&t(z) = z where TeM*. Then 

(5.16) w'(I) is invertible <*=^> kcr(A/ 2 (T) - Id) = Span (JVFi(z), . . . , JVF n (z)). 

(5.17) < w'{I)- 1 (w(I)), w(I) <=> JVH(z) i (M Z (T) - Id)(T z E E ). 

Suppose moreover |5.7| ) and \5.10\) . Then periodic orbits of Tie x suppf are Y^E-normal. Besides, 
if (T, z) is such a periodic point, then we have E\(T, z) = ker(M z (T) — Id) 2 = M. 2n . 

Note that Lemma 1531 actually follows from Proposition 14.51 

Proof. Let us denote by 

(5.18) W(z) :=Span (VF^z), . . . ,VF n (z)). 

For any integrable system of the form (|5.ip . we have Ei(T,z) = R 2n . To see this, just apply 
Lemma l2~5l with V = W = W(z). We introduce the canonical projection 

(5.19) u : R" -> IT. 

Of course, the differential of u at any point of R" is an isomorphism. Note that in (|5.4[) . to 
be rigorous, we should say 9t(9,I) = dgu(V H(I)), where 9 £ l n is such that u(9) = 9. If 
(9,1) := ip(z), we denote by M {9J) (t) := (fi' T (9,I). Then for all (t,x) E TgT" x K" we have 

(5.20) (M iej) (T) - Id)(r,x) =T(d § u(H"(I)x),0). 

We can read on this last formula that TgT 1 x {0} C ker(M (e /) (T) - Id)) and that H"(I) is 
invertible if and only if dim(kcr(A/(g j)(T) — Id)) = n. 
Using properties of ijj, we note that 

(5.21) = T n x £ E , d z 4>(JVH(z)) = (d g u(VH(I)),0). 

Thus, in view of (pT2T)j) . JVH(z) G (M Z (T) - Id)(T z E B ) if and only if there exists (t,x) in 
T e T n x (M.VH(I)) 1 - such that (d § u(VH(I)), 0) = (M (ej) (T) - Id)(r,x) = T(d § u(H" (I)x), 0), 

i.e. ifV-ff(I) = TH"(I)x, where x ±VH(I). That exactly means < ^(i")" 1 («;(/)), >= 0. 
Since M Z (T) is symplectic we have 

(5.22) Im (M Z (T) — J) = ker(* M Z (T) - /d) 1 " = J[ker(M z (T) - Id)]- 1 . 

We always have JW(z) C ker(Af 2 (T) — Id). Now if we make assumption (|5.7p . in view of 
(I5.10|) . we have equality. Moreover due to the involutions, we have W(z) C JW(z) ± , and in fact 
W(z) = JW(z) 1 - for dimensional reasons. Coming back to (|5.22p . we obtain Im (M Z (T) — I) = 
JW(z) = ker(M z (T) - Id), which means precisely (M Z (T) - Id) 2 = 0. □ 

We can now give a rigorous mathematical proof of the Berry- Tabor trace formula (see original 
papers [5] and [H]). 
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Theorem 5.4. Make hypotheses H5.1]) . (5. Sty , H5.S]) . (5. 7\ ) and i5.10]) . Suppose for clarity that 
^ suppf. Then 

M:= T ^ n gZ" 

where [3i t T is an integer such that e l ^^ I - T = i n ~ x sign(K (7)) , is i/ie Gaussian curvature of 

Tie at I £ R" (cf \5.9\) and (5.11]) ), and w(I) is the frequency function on the torus I (cf (5.6]) ). 



Proof. In view of Lemma 15.31 and Proposition 14. 5[ we can apply Theorem 12.41 Here, according 
to section 15.11 connected components of G E are one to one with the set of (T, I) e supp/ x 
such that Tw(I) € (27rZ) ra . Each connected component has the dimension of a periodic torus, 
i.e. equal n. We obtain modulo 0(h°°): 
(5.23) 



9 E {h) = (2ith)^^(E) 



E W )/(T)i- I d(T, z)da 1 {z) + £ /,'„,. 

E K Xsunpf: \ 



(7,T)6S B Xsupp/: 

where do - / is the Riemannian measure on the torus Ti := ^ -1 ({(0, /)}) and A(T,I) is the 
constant value of the action on {T} x L/. 

First of all, we simplify the formula of the density. Let (T, z) be a periodic point of supp/ xSg. 
Since E x = R 2 ", we have det(M z (T) - Id)\ Vi =1 and det(w | Bl ) = 1. We introduce 

£ 4 (T,z) -g^T,*)^ 

For sake of simplicity in the following ,we omit '(T, z)' in '£j(T, z)' notation. First we note that 
we have £4 = J£i = J£2 = W(z) (with the notation of (|5.18p ). We claim that 

J(M Z (T) - Id)8. 4 = £4. 

Indeed we have seen in the proof of Lemma [O] that (M Z (T) - Id)(M. 2n ) = JW(z) and R 2 ™ = 
£1 + £4 with £1 = ker(Af z (T) — Id). We note that since £5 = £4 n T z Y,e, we have 

£ 5 ©RViJ(z) = £4. 

Let £1 € £4 such that (M Z (T) — Id)ei = JVH(z). Let j3" be an orthonormal basis of £5. Then, 
since JVH(z) £ (M Z (T) - Id){T z T, E ), P ■= (ei,/3") is a basis of £ 4 . Let f3' := (VH(z),p"). It 
is of course also a basis of £4. Then 

det[J{M z (T) - Id) hi ] = det[J(M z (T) - Id)| e J/9,/J' det[/d]^, . 
Since V-Ef(z) _L /3", we have det[J(M z (T) - /d)| £ J/3,/3' = - det(n £5 J(M Z (T) - Id)\ ) and 
det[Id]p> i0 = \\VH(z)\\ 2 / < e u VH(z) >. Thus 

\\VH(z)\\ 2 det(n £s J(A/ 2 (T) -/d)k 8 ) = - < ei,Vif(«) > dct[(M z (T) - Jd) fi J. 
Therefore, in view of Theorem 12.41 

(5-24) d(T, z) 2 = ; ; K -r4 ^ ; r 

K ' 1 ' <e 1 (T,z),VH(z)>det(J(M z (T)-I) k J 

where ei(T, z) is the element of £4 such that (M z {T)-Id)[ei(T, z)] = JVH(z). Let a = (0, 1) := 
ip(z) and M a (T) := d^M^T^d^' 1 . As £4 = J£i, and as d z i/)(£i) = T g T n x {0}, we have 

det(J(M z (T) - I) h J - dct((M z (T) - I)J Wi ) = det[(M Q (T) - Id^Jd^- 1 )^^]. 

Using notations of ([519]) we introduce G e : {0} x R" -> T g T n x {0} defined by 

G e {0,x) := (d $ u(x),0). 
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We have 

(5.25) det(J(M z (T) - = dct[G^(M a (T) - Id^^dctKd^Jd^Gg)^^}. 
We introduce i : {0} x R" — > R" defined by i(0, x) := x. Of course, 

det[Gf(M a (T) - /d) Uo}XB „] = det [iG~\M a {T) - Id^ 1 }. 
In view of lpT20j) . for all iG^ l (M a {T) - Id)^ 1 = TH"(I). Therefore we have 

(5.26) det[Ge\M a (T)-Id) ko}XKn ]=T n det(H"(I)). 

We obtain a paramctrisation of L/ = ?/> -1 ((0, /)) with the following fonction ip : [0, 2n] n — > L/ 
defined by 

(5.27) i-^Mt),/), Vre[0,2^]". 
Differentiating the last equation leads to 

d^tp o d T <£ oi = G e . 

Thus 
(5.28) 

det[(d l/3(r) '0oJod I p( T )V' _lo Ge)| {o}xR „] = det[(d ¥ , (T) ^o Jo^^oi)^^] = det[(iod v(T )V>oJ<£ r ¥>)| Rn ]. 

We claim that we have i o d v ( T \ip =* d r y> o J on £4. To check this, let /? G £4 and y G K™. By 
definition of t d T ip, we have 

<* d T ipoJf3,y >=< Jf3,d T ip(y) >= w lp{v{T)) {d ip{T) -ip(p),d v{T] 'ip o d T ip{y)). 

where w) is the symplectic form on T" x R n . We have used here that ip is a symplectic map. 
Now, in view of (|5.27|) . and using the definition of w, the last term of the preceeding equation is 
equal to 

wVMt))(( c W)'0(/3), (d T u(y),0)) =< io d v{r) ip(f3),y > . 
Thus we have proven that i o d^^ip =* d T (p o J on £4. Coming back to ()5.28p . we obtain 

det[(d v(T) V>o J o d^i/)- 1 o G e )\ {0}XKn } = (-1)™ det(< 'd T <pd T cp) =: (-l) n g v ( T ). 

In view of (|5.25[) and (|5.26p . we have, if z = <p(t) 

(5.29) det(J(M z (T) - J)| £ J = T n det( J ff"(/))^(r)(-l) n . 
We claim that 

(5.30) < ei,VfT(z) >= i < H"{I)- l VH(I),VH{I) > . 
Indeed by definition of £1, we have in view of (|5.21[) 

(M Z (T) - Id)d z tP(e 1 ) = (d § u(VH(I)),0) 

Let (ti,xi) := d z ^{e x ). Using ([5^0]) we obtain Tx x = H"(I)~ 1 VH(I). Now, -0 being a 
symplectic map, we have 

<e u VH{z) >=< Je lt JVH(z) >= w^ z) (d t i/)(Ei), (dg«(V#(J)), 0)) =< a*, Vlr(J) >, 

which proves (|5.30p . In view of (|5.24[) . (|5.29|) and (|5.30[) . we have obtained if f{r) = z 

_7>,--(n+i) 1 

d(T,z) 2 = ^ 7v 

< H"{I)- l VH{I),VH{I) > T n det{H"{I)) 9 V ( T ) 

In view of Lemma 15.21 we obtain 



(_1 y i ? '-( n +i) 1 
\w(I)\ 2 K(I)\M\ n 1 ^(r) 
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where 2ttM := Tw(I). There exists an integrer (i = (3{T,I) such that = H n+1 )sign(if (I))(-1) T 
and that if ip(r) = z, then 



d(T,z). g |(r) = (2 7 r)-("- 1 ). 



e if/9 



w(I)\^\K(I)\\M\^- 

1 f If i e iJ i p 

d(T,z)dai(z) = — d(T,(p{T))g${T)dT 



Note that, since z i— ► d(T, z) is continuous, our integer (3 is constant on the connected torus T/. 
Endly, computing in action-angles coordinates gives 

A(T, I) = 2n <M;I >= T < w(I);I > . 

This ends the proof of Theorem 15.41 □ 
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